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'. Abstract 
O ■ 

' A detailed quantitative analysis of the system (ppe) placed in magnetic field ranging from — 

(N : 

Q ■ 4.414 X 10^'^ G is presented. The present study is focused on the question of the existence of the 

Q . molecular ion H2 in a magnetic field. As a tool, a variational method with an optimization of 

I the form of the vector potential (optimal gauge fixing) is used. It is shown that in the domain 

of applicability of the non-relativistic approximation the system (ppe) in the Born-Oppenheimer 

• approximation has a well-pronounced minimum in the total energy at a finite interproton distance 

, for B < 10^^ G, thus manifesting the existence of H2 ■ For B > 10^^ G and large inclinations 

(N ■ ^ 

T— I ' (of the molecular axis with respect to the magnetic line) the minimum disappears and hence the 

(N : 

■ molecular ion H2 does not exist. It is shown that the most stable configuration of H2 always 

^ ■ + 

Q^. corresponds to protons situated along the magnetic line. With magnetic field growth the ion H2 

O I becomes more and more tightly bound and compact, and the electronic distribution evolves from 

^ ' a two-peak to a one-peak pattern. The domain of inclinations where the H2 ion exists reduces 

^ ■ with magnetic field increase and finally becomes 0° — 25° at B = 4.414 x 10^'^ G. Phase transition 

type behavior of variational parameters for some interproton distances related to the beginning of 
the chemical reaction H2 H + p is found. 
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I. INTRODUCTION 



Many years have passed since the moment when theoretical quahtative arguments were 
given that show that in the presence of a strong magnetic field the physics of atoms and 
molecules exhibits a wealth of new, unexpected phenomena even for the simplest systems 
|l|, 1^. In particular, a chance that unusual chemical compounds may be formed which do not 
exist without magnetic field was mentioned. In practice, the atmosphere of neutron stars, 
which is characterized by the presence of enormous magnetic fields 10^^ — 10^'^ G, as well 
as other astronomical objects carrying large magnetic fields (> 10^ G) provide a valuable 
paradigm where this physics could be realized. Recently, the experimental data collected by 
the Chandra X-ray observatory revealed certain irregularities in the spectrum of an isolated 
neutron star 1E1207.4-5209. These irregularities can be interpreted as absorption features 
at ~ 0.7 KeV and ~ 1.4 KeV of possible atomic or molecular nature 0. 

One of the first general features observed in standard atomic and molecular systems placed 
in a strong magnetic field is an increase of both total and binding energies, accompanied by a 
drastic shrinking of the electron localization length in both the longitudinal and transverse 
directions. Naturally, this leads to a decrease of the equilibrium distance with magnetic 
field growth. This behavior can be considered to be a consequence of the fact that for large 
magnetic fields the electron cloud takes a needle-like form extended along the magnetic field 
direction and the system becomes effectively quasi-one- dimensional It is obvious that 
the phenomenon of quasi-one-dimensionality enhances the stability of standard atomic and 
molecular systems from the electrostatic point of view. In particular, molecules become 
elongated along the magnetic line forming a type of linear molecular polymer (for details 
see the review articles @, H). It also hints at the occurrence of exotic atomic and molecular 
systems which do not exist in the absence of a magnetic field. Motivated by these simple 
observations it was shown in Refs. |P, |^ that three and even four protons can be bound by 
one electron. This shows that exotic one-electron molecular systems H^'^ and H^'^ can exist 
in sufficiently strong magnetic fields in the form of linear polymers. However, the situation 
becomes much less clear (and also much less studied) when the nuclei are not aligned with 
the magnetic field direction, and thus in general do not form a linear system. Obviously, 
such a study would be important for understanding the kinetics of a gas of molecules in 
the presence of a strong magnetic field. As a first step towards such a study, even the 
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simplest molecules in different spatial configurations deserve attention. Recently, a certain 
spatial configuration of H^'^ was studied in detail P). It was shown that in the range of 
magnetic fields 10^ < B < 10^^ G the system (pppe), with the protons forming an equilateral 
triangle perpendicular to the magnetic lines, has a well-pronounced minimum in the total 
energy for a certain size of triangle. The goal of the present work is to attempt for the 
first time to carry out an extensive quantitative investigation of the ground state of H2 in 
the framework of a single approach in its entire complexity: a wide range of magnetic field 
strengths (0 — 4.414 x 10^^ G), arbitrary (but fixed) orientation of the molecular axis with 
respect to the magnetic line and arbitrary internuclear distances. We are going to carry out 
this study in the Born-Oppenheimer approximation at zero order - assuming protons to be 
infinitely heavy charged centers. In principle, when the molecular axis is perpendicular to 
the magnetic line the system (ppe) acquires extra stability from the electrostatic point of 
view. Electrostatic repulsion of the classical protons is compensated for by the Lorentz force 
acting on them. 

It is well known that the molecular ion H2 is the most stable one-electron molecular 
system in the absence of a magnetic field. It remains so in the presence of a constant 
magnetic field unless B > 10^^ G, where the exotic ion appears to be the most bound 
(see [0]). The ion H2 has been widely studied, both with and without the presence of a 
magnetic field, due to its importance in astrophysics, atomic and molecular physics, solid 



state and plasma physics (see QI-ISC] and references therein). The majority of the previous 
studies were focused on the case of the parallel configuration, where the angle between 
the molecular axis and the magnetic field direction is zero, 6 = 0°. The only exception 



is Ref. |]T6[, where a detailed quantitative analysis was performed for any 9 but for a single 
magnetic field B = 1 a.u. . Previous studies were based on various numerical techniques, but 
the overwhelming majority used different versions of the variational method, including the 
Thomas- Fermi approach. As a rule, in these studies the nuclear motion was separated from 
the electronic motion using the Born-Oppenheimer approximation at zero order (see above). 
It was observed at the quantitative level that magnetic field growth is always accompanied 
by an increase in the total and binding energies, as well as a shrinking of the equilibrium 
distance. As a consequence it led to a striking conclusion about the drastic increase in the 
probability of nuclear fusion for H2 in the presence of a strong magnetic field [Q. 

In the present study we will also use a variational method. Our consideration will be 
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limited by a study of the Ip-state, which reahzes the ground state of the system if the bound 
state exists [0]. We will construct state-of-the-art, non-straightforward, 'adequate' trial 
functions consistent with a variationally optimized choice of vector potential. We should 
stress that a proper choice of the form of the vector potential is one of the crucial points 
which guarantee the adequacy and reliability of the consideration. In particular, a proper 
position of the gauge origin, where the vector potential vanishes, is drastically important, 
especially for large interproton distances. For the parallel configuration, 6 = 0° the present 
work can be considered as an extension (and also an development) of our previous work 



2^ . It is necessary to emphasize that we encounter several new physical phenomena which 



occur when the molecular axis deviates from the magnetic field direction. If the magnetic 
field is sufficiently strong, B > 10^^ G and inclination 6 is larger than a certain critical 



angle the ion H2 does not exist in the contrary to a prediction in Refs.|ll], |p, |24|. This 



prediction was based on an improper gauge dependence of the trial functions, which caused 
a significant loss of accuracy and finally led to a qualitatively incorrect result. We find that 
in the weak field regime the [ppe) system in the equilibrium position at any inclination, 
the electronic distribution peaks at the positions of the protons. While at large magnetic 
fields the electronic distribution is characterized by single peak at the midpoint between 
two protons. This change from a two-peak to a one-peak configuration appears around 
B ~ 10^'^ — 10^^ G with a slight dependence on the inclination angle 6. From a physical 
point of view the former means that the electron prefers to stay in the vicinity of a proton. 
This can be interpreted as a dominance of the if-atom plus proton configuration. The latter 
situation implies that the electron is 'shared' by both protons and hence such a separation 
to iY-atom plus proton is irrelevant. Therefore, we can call the two-peak situation 'ionic' 
coupling, while the one-peak case can be designated as 'covalent' coupling, although this 
definition differs from that widely accepted in textbooks (see, for example |31|). Thus, we 
can conclude that a new phenomenon appears - as the magnetic field grows the type of 
coupling changes from 'ionic' to 'covalent'. At large internuclear distances the electron is 
always attached to one of the charged centers, so the coupling is 'ionic'. 

One particular goal of our study is to investigate a process of dissociation of the {ppe) 
system: H2 H + p which appears with increase of interproton distance. It is clear from a 
physical point of view that at large distances the electronic distribution should be first of the 
two-peak type and then should change at asymptotically large distances, to a single-peak 
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one, but with a peak at the position of one of the protons. Somehow this process breaks 
permutation symmetry and we are not aware of any attempt to describe it. In our analysis 
this phenomenon appears as a consequence of a change of position of the gauge origin with 
increase of interproton distance. 

From the physical point of view it is quite interesting to note how the (ppe) system 
behaves at very large interproton distances. This domain is modelled by an iJ-atom plus 
proton interaction. The interaction corresponds to (magnetic-field-inspired-quadrupole) + 
charge interaction and is dominant comparing to the standard Van der Waals force. For 
small inclinations the above interaction is attractive as in the Van der Waals case, but 
becomes repulsive for large inclinations. This implies that the potential curves approach to 
asymptotic value of total energy at the large interproton distances from above contrary to 
the Van der Waals case. 

The Hamiltonian which describes two infinitely heavy protons and one electron placed in 
a uniform constant magnetic field directed along the 2;— axis, B = (0, 0, B) is given by (see 
e.g. 
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n = p^ + ^ - - - - + {pA + Ap) + A\ (1) 
R Ti r2 

(see Fig.l for notations), where p = —iV is the momentum, ^ is a vector potential, which 
corresponds to the magnetic field B. Hence, the total energy Et of H2 is defined as the 
total electronic energy plus the Coulomb energy of proton repulsion. The binding energy is 
defined as an affinity to have the electron at infinity, = B — Et- The dissociation energy 
is defined as affinity to have a proton at infinity, Ed = Eh — Et-, where Eh is the total 
energy of the hydrogen atom in magnetic field B. 

Atomic units are used throughout {h=me=e=l) albeit energies are expressed in Rydbergs 
(Ry). Sometimes, the magnetic field B is given in a.u. with Bq = 2.35 x 10^ G 



II. OPTIMIZATION OF VECTOR POTENTIAL 



It is well known that the vector potential for a given magnetic field, even in the Coulomb 
gauge (V ■ ^) =0, is defined ambiguously, up to a gradient of an arbitrary function. This 
gives rise a feature of gauge invariance: the Hermitian Hamiltonian is gauge-covariant, while 
the eigenenergies and other observables are gauge- independent. However, since we are going 
to use an approximate method for solving the Schroedinger equation with the Hamiltonian 
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FIG. 1: Geometrical setting for the H2 ion placed in a magnetic field directed along the z-axis. 
The protons are situated in the y — z plane at a distance R from each other and marked by bullets. 
O is the origin of coordinates which is chosen to be on the bold-dashed line which connects protons; 
O'(0, y, Z) is the mid-point between protons. It is assumed that the gauge center coincides with 
O. 00' measures the distance between the gauge center and the mid-point between the proton 
positions (see text and eq.(4)). 

(0), our approximation of eigenenergies can well be gauge- dependent (only the exact ones 
are gauge-independent). Hence one can choose the form of the vector potential in a certain 
optimal way. In particular, if the variational method is used, the vector potential can be 
considered as a variational function and be chosen by a procedure of minimization. 

Let us consider a certain one-parameter family of vector potentials corresponding to a 
constant magnetic field B = (0, 0, B) 



where ^ is parameter, in the Coulomb gauge. The position of the gauge center or gauge 
origin, where A{x, y, z) = 0, is defined hj x = y = 0, with z arbitrary. For simplicity we fix 
2; = 0. If^ = l/2we get the well-known and widely used gauge which is called symmetric or 



A = Bi{^-l)y, ^x, 0) , 



(2) 
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circular. If ^ = or 1, we get the asymmetric or Landau gauge (see 0]). By substituting 
(0) into (|I|) we arrive at a Hamiltonian of the form 



n = -V' + ^ ~ - - - - 2tB[{^ - l)yd, + ^xdy] + B^iex^ + (1 - OV] , (3) 
n ri r2 

where R is the interproton distance (see Fig. 1). 

It is evident that for small interproton distances, R, the electron prefers to be near the 
mid-point between the two protons (coherent interaction with the protons). In the opposite 
limit, R large, the electron is situated near one of the protons (this is an incoherent situation 
- the electron selects and then interacts essentially with one proton). This fact, together 
with naive symmetry arguments, leads us to a natural assumption that the gauge center is 
situated on a line connecting the protons. Therefore, the coordinates of mid-point between 
protons are 

Y = ^sme , Z = ^cose , (4) 

(see Fig.l), where is a parameter. Thus, the position of the gauge center is effectively 
measured by the parameter d - a relative distance between the middle of the line connecting 
the protons and the gauge center. If the mid-point coincides with the gauge center, then 
d = 0. On other hand, if the position of a proton coincides with the gauge center, then d = 1 
or d = —1. Hence the parameter d makes sense as a parameter characterizing a gauge. 

The idea of choosing an optimal (convenient) gauge has been widely exploited in quantum 
field theory calculations. It has also been discussed in quantum mechanics and, in particular, 
in connection to the present problem. Perhaps, the first constructive (and remarkable) 
attempt to realize the idea of an optimal gauge was made in the eighties by Larsen |TT| . In 
his variational study of the ground state of the H2 molecular ion it was explicitly shown that 
for a given fixed trial function the gauge dependence of the energy can be quite significant. 
Furthermore, even an oversimplified optimization procedure improves the accuracy of the 



numerical results |^| 



Our present aim is to study the ground state of (d) or, more concretely, (H). We propose 
a different way of optimization of vector potential than those discussed by previous authors. 
It can be easily demonstrated that for a one-electron system there always exists a certain 
gauge for which the ground state eigenfunction is a real function. Let us fix a vector 
potential in (|l|). Assume that we have solved the spectral problem exactly and have found 



the ground state eigenfunction. In general, it is a certain comj^/ex function with a non-trivial, 
coordinate-dependent phase. Treating this phase as a gauge phase and then gauging it away, 
finally results in a new vector potential. This vector potential has the property we want 
- the ground state eigenfunction of the Hamiltonian (0) is real. It is obvious that similar 
considerations are valid for any excited state. In general, for a given eigenstate there exists 
a certain gauge in which the eigenfunction is real. For different eigenstates these gauges 
can be different. It is obvious that similar situation occurs for any one-electron system in a 
magnetic field. 

Dealing with real trial functions has an immediate advantage: the expectation value of 
the terms proportional to A in (|l]) (or ~ 5 in (^)) vanishes when it is taken over any real, 
normalizable function. Thus, without loss of generality, the term ~ in (0) can be omitted. 
Thus, we can use real trial functions with explicit dependence on the gauge parameters ^ 
and d. These parameters are fixed by performing a variational optimization of the energy. 
Therefore, as a result of the minimization we find both a variational energy and a gauge for 
which the ground state eigenfunction is real, as well as the corresponding Hamiltonian. One 
can easily show that for a system possessing axial (rotational) symmetry |^| the optimal 



gauge is the symmetric gauge ^ = 1/2 with arbitrary d. This is precisely the gauge which 
has been overwhelmingly used (without any explanations) in the majority of the previous 
research on in the parallel configuration [|T| - [^. However, this is not the case if 7^ 0°. 
For the symmetric gauge the exact eigenfunction now becomes complex, therefore complex 
trial functions must be used. But following the recipe proposed above we can avoid complex 
trial functions by adjusting the gauge in such a way the eigenfunction remains real. This 
justifies the use of real trial functions. Our results (see Section IV) lead to the conclusion 
that for the ground state the optimal gauge parameter varies in the interval ^ G [0.5 , 1]. 

III. CHOOSING TRIAL FUNCTIONS 

The choice of trial functions contains two important ingredients: (i) a search for the gauge 
leading to the real, exact ground state eigenfunction and (ii) performance of a variational 
calculation based on real thai functions. The main assumption is that a gauge corresponding 
to a real, exact ground state eigenfunction is of the type @ (or somehow is close to it) 



In other words, one can say that we are looking for a gauge of type (p|) which admits the 
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best possible approximation of the ground state eigenfunction by real functions. Finally, in 
regard to our problem, the following recipe of variational study is used: As the first step, 
we construct an adequate variational real trial function \E'o ]^/ , for which the potential 
Vq = reproduces the original potential near Coulomb singularities and at large distances, 
where ^ and d would appear as parameters. The trial function should support symmetries of 
an original problem. We then perform a minimization of the energy functional by treating 
the free parameters of the trial function and C,,d on the same footing. In particular, such an 
approach enables us to find the optimal form of the Hamiltonian as a function of ^, d. 

The Hamiltonian (|l|) gives rise to different symmetry properties depending on the ori- 
entation of the magnetic field with respect to the internuclear axis. The most symmetric 
situation corresponds to 6 = 0°, where invariance under permutation of the (identical) 
charged centers P : (1 ^ 2) together with : {z —z) holds. Since the angular momen- 
tum projection £^ = m is conserved, accounts also for the degeneracy m —m. Thus, we 
classify the states as l<Jg,u, ■ ■ ■ l7i"g,M, '^'^9,u ■ ■ ■ ^Sg,u, '^5g,u ■ ■ ■, where the numbers 1,2,... 
refer to the electronic states in increasing order of energy. The labels a, vr, 5 . . . are used 
to denote |m| = 0,1,2..., respectively, the label g (u) gerade (ungerade) is assigned to 
the states of even (odd) parity P of the system. At 6 = 90" the Hamiltonian still remains 
invariant under the parity operations P and Pz, while the angular momentum projection is 
no longer conserved and m is no longer a quantum number. The classification in this case 
is 1^„, 2^^, . . ., where the sign +(— ) is used to denote even (odd) 2;-parity. Eventually, for 
arbitrary orientation, only parity under permutations P is conserved. In general, we refer 
to the lowest gerade and ungerade states in our study as Ig and 1„. This is the only unified 
notation which make sense for all orientations 0° < < 90°. 

The above recipe (for the symmetric gauge where ^ = 1/2, = 0) was successfully apphed 
in a study of the if2^-ion in a magnetic field for the parallel configuration 6 = 0" and also 
for general one-electron linear systems aligned along the magnetic field [0. In particular, 
this led to the prediction of the existence of the exotic ions at B > 10^° G and m a 
linear configuration at B > 10^'^ G 0, |^. Recently, this recipe was used for the first 
time to make a detailed study of the spatial configuration It was demonstrated 

that inconsistency between the form of vector potential and a choice of trial functions can 
lead to non-trivial artifacts like existence of spurious bound states (see ||2(j||). 

One of the simplest trial functions for Ig state which meets the requirements of our 
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criterion of adequacy is 



(cf. p2| , p^), where ai, /^i^,. and (3iy are variational parameters and ^ is the parameter 
of the gauge (2). The first factor in the function (|^), being symmetric under interchange 
of the charge centers ri r2, corresponds to the product of two Is-Coulomb orbitals 
centered on each proton. It is nothing but the celebrated Heitler-London approximation 
for the ground state lag. The second factor is the lowest Landau orbital corresponding 
to the vector potential of the form Eq. (||). So, the function (|^) can be considered as a 
modification of the free field Heitler-London function. Following the experience gained in 
studies of H2 without a magnetic field it is natural to assume that Eq. (j^) is adequate 
to describe interproton distances near equilibrium. This assumption will be checked (and 
eventually confirmed) a posteriori, after making concrete calculations (see Section IV). 
The function (13) is an exact eigenfunction in the potential 



Vi = ^ = 2al - 2B[(3,,^ + A,(l - 0] + ^B'[(3iex' + f3l{l - OV] + 2«?(ni ■ ^2 



+4q!i5 



- 2m 



1 1 

ri r2 



ri r2 
where yi^2 are the y-coordinates of protons (see Fig.l). The potential Vi reproduces the 
functional behavior of the original potential (^ near Coulombic singularities and at large 
distances. These singularities are reproduced exactly when = Piy = 1/2 and ai = 1. 

One can construct another trial function which meets the requirements of our criterion 
of adequacy as well, 



(cf. p2|, |2^). This is the celebrated Hund-Mulliken function of the free field case multiplied 
by the lowest Landau orbital, where 02, f32x and (32y are variational parameters. From 
a physical point of view this function has to describe the interaction between a hydrogen 
atom and a proton (charge center), and, in particular, models the possible dissociation mode 
of H2 into a hydrogen atom plus proton. Thus, one can naturally expect that for sufficiently 
large internuclear distances R this function prevails, giving a dominant contribution. Again 
this assumption will be checked a posteriori, by concrete calculations (see Section IV). 
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There are two natural ways to incorporate the behavior of the system in both regimes - 
near equihbrium and at large distances - into a single trial function. It is to make a linear 
or a nonlinear interpolation. The linear interpolation is given by a linear superposition 

^^3^ = v4i^i + A2^2 , (7) 

where Ai or A2 are parameters and one of them is kept fixed by a normalization condition. 
In turn, the simplest nonlinear interpolation is of the form 



(cf. IP^ , where as, ^4, /S^x and /S^y are variational parameters. This is a Guillemin- 
Zener function for the free field case multiplied by the lowest Landau orbital. If 0:3 = q;4, 
the function (|) coincides with (|). If 0:4 = 0, the function (H) coincides with (|). 
The most general Ansatz is a linear superposition of the trial functions (0) and (P), 

^ = Ai^i + A2^2 + ^3^36 , (9) 

where we fix one of the A's and let all the other parameters vary. Finally, the total number 
of variational parameters in (^), including R,C,,d, is fifteen for the ground state. For the 
parallel configuration, 6 = 0°, the parameters C, = 0.5, d = are fixed in advance and also 
Pix = Piy,P2x = P2y,P3x = Psy Heucc the uumbcr of free parameters is reduced to ten for 
the ground state. Finally, with the function (^) we intend to describe the ground state for 
all magnetic fields where non-relativistic consideration is valid, B < 4.414 x 10^'^ G, and for 
all orientations of the molecular axis. 

Calculations were performed using the minimization package MINUIT from CERN-LIB. 
Numerical integrations were carried out with a relative accuracy of ~ 10^^ by use of the 
adaptive NAG-LIB (DOIFCF) routine. All calculations were performed on a PC Pentium-Ill 
800 MHz. 



IV. RESULTS 

We carry out a variational study of the system (ppe) with infinitely heavy protons in 
the range of magnetic fields < i? < 4.414 x 10^'^ G, inclinations 0° — 90°, for a wide 
range of interproton distances. For magnetic fields B < 10^^ G the system displays a well- 
pronounced minimum in the total energy at all inclinations. However, for B > 10^^ G at 
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large inclinations the minimum in the total energy disappears, while for small inclinations a 
minimum continues to exist. This picture describes the domain of existence of the molecular 
ion H2 ■ In general, we confirm a qualitative results by Khersonskij about the non- 
existence of a minimum at finite distances on the total energy surfaces of the system (ppe) 
at sufficiently strong magnetic fields for some far from parallel orientations. It is worth 
mentioning that the variational study in |Q was carried out with a trial function somewhat 
similar to that of Eq.(|^), which, however, does not fully fulfill our criterion of adequacy. 
The potential corresponding to this function reproduces correctly the original potential near 
Coulomb singularities and ~ growth at large distances. However, it generates growing 
terms ~ p which implies a reduction of the rate of convergence of a perturbation theory for 
which the variational energy represents the first two terms (see the discussion in [^). Also, 
this trial function is not satisfactory from the point of view of gauge invariance. However, 
in spite of all the above-mentioned deficiencies it led to qualitatively correct picture. 

In Figs. 2-5 the total energy Et of the {ppe) system as a function of interproton distance 
R for several values of the magnetic field strength and different values of inclination 9 is 
shown. For magnetic fields B < 10^^ G and for all inclinations 0° — 90°, each plot displays 
a well-pronounced minimum at i? = -Rgg, manifesting the existence of the molecular system 
H2 ■ For B = la.u. and R < 3.5 a.u. (see Fig. 2) our results are similar to the results 
of pi| , P^ ] - for fixed R the potential energy Et grows with inclination. In general, at 
large R > Reg and for 6 > 0° all the curves behave alike: they reveal a maximum and 
then tend (from above) to the total energy of the hydrogen atom. For 6 = 0° the potential 
curves approach to the asymptotics from below, displaying in general a behavior similar 
to the field-free case, to a Van der Waals-force-inspired behavior. This behavior is related 
to the fact that at large R the configuration if-atom + proton appears. The if-atom has 
quadrupole moment, Q B^ (see p^, 0). Hence at large distances the total energy 



is defined by a quadrupole moment - charge interaction 

eQiB)P2icose) 
Et = , (10) 

where P2 is the second Legendre polynomial. At small inclinations P2{cos6) is positive, the 
total energy is negative, thus corresponding to attraction between the quadrupole and the 
charge. Therefore, the total energy curve approaches to the asymptotics from below. For 
large inclinations P2(cos 9) is negative and the total energy is positive. Thus, this corresponds 
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FIG. 2: Total energy Et in Rydbergs of the (p|?e)-system as function of interproton distance R for 
different inclinations at B = 2.3505 x 10^ G ( la.u.) 



to repulsion between quadrupole and charge, and implies an existence of maximum of the 
total energy for large interproton distances R > Req. We observe the maximum in all 



Figs. 2-5. It is worth mentioning that in the calculations 0, for B = la.u. and 6 = 90° 
(and other inclinations) the maximum was not observed in contradiction to our predictions 
(see Fig. 2 and also below Fig. 9). Looking at Fig. 2 it is interesting to compare a rate 
with which potential curves are approaching to the asymptotic total energy at large R. This 
asymptotic energy is equal to the total energy of the hydrogen atom, Eh = —0.6623Ry, while 
E^='^°{R = 8a.u.) = -0.6647 (from below), E^=^^'' {R = 8a.u.) = -0.6576 (from above), 
E^^^°{R = 8 a.u.) = —0.6620 (from above). Thus, any deviation does not exceed 1%. There 
exists a different manner of viewing these results. It can be treated as a demonstration of 
the quality of our trial function (|^) but for the calculation of the total energy of the atom 

(0- 

However, the situation is drastically different for B > 10^^ G, see Figs. 4-5. There exists 
a certain critical angle 6cr, such that for 6 < 9cr the situation remains similar to that given 
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FIG. 3: Total energy Et in Rydbergs of the (ppe)-system as function of interproton distance R for 
different inclinations at B = 10^^ G 



above - each potential curve is characterized by a well-pronounced minimum at finite R. 
With increase of the inclination, at > 6cr the minimum in the total energy first becomes 
very shallow with Et > Eh and ceases to exist at all. We were unable to localize with 
a confidence the domain in R which correspond to a shallow minimum which leads to the 
possible dissociation H2 H + p that was predicted in [|ll|] as well as being discussed in 



our previous work pj] . We consider that the prediction of dissociation for large inclinations 
emerged as an artifact of an improper choice of the gauge fixing (see the discussion above). 
A detailed study of the transition domain (existence non-existence) of H2 is beyond 
of scope of the present article. In any case, such a study requires much more accurate 
quantitative techniques as well as a sophisticated qualitative analysis. Schematically, the 
situation is illustrated in Fig. 6. 

It is quite interesting to explore the variation of the vector potential (^ for 6 7^ 0°, in 
particular, the position of the gauge center as a function of interproton distance R and 



magnetic field p9[. In Figs. 7 a,b for B = la.u. and Figs. 8 a,b for 10 G, correspondingly, 
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FIG. 4: Total energy Et in Rydbergs of the (ppe)-system as function of interproton distance R for 
different inclinations at i? = lO"*^^ G. The result by Wille [14| is shown by bullet (see text) 



both the ^— and d— dependence are presented (see (H) and discussion in Sec. III). This 
dependence is very similar for all magnetic fields studied. It is worthwhile to emphasize that 
for all the potential curves given the minimum (in other words, the equilibrium position) 
at -R = Reg somehow corresponds to a gauge close to the symmetric gauge: ^ > 1/2 |^ 
and Y = Z = 0{d = 0). A similar situation holds for small interproton distances, R < Req. 
However, for large R,R> Reg the parameter ^ grows smoothly, reaching a maximum near 
the maximum of the potential curve which we denote hj R = Rcr- It then falls sharply 
to the value ~ 1/2. In turn, the parameter d remains equal up to i? = Rcr (which 
means the gauge center coincides to the mid-point between protons), then sharply jumps to 
1 (gauge center coincides with the position of a proton), displaying a behavior similar to a 
phase transition. It is indeed a type of phase transition behavior stemming from symmetry 
breaking: from the domain R < Rcr, where the permutation symmetry of the protons holds 
and where the protons are indistinguishable, to the domain R > Rcr, where this symmetry 
does not exist and electron is attached to one particular proton. Such a type of 'phase 
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FIG. 5: Total energy Et in Rydbergs of the (ppe)-system as a function of interproton distance 
R for different inclinations at = 10^^ G. Plots for = 45°, 60", 90° consist of two parts - (i) 
solid line, when d is kept fixed, d = (gauge center coincides with the mid-point between protons) 
and which displays a minimum, and (ii) the dotted line is the result of minimization, when the 
parameter d is released. 

transitions' is typical in chemistry and is called a 'chemical reaction'. Hence the parameter 
Rcr characterizes a distance at which the chemical reaction H2 ^ H + p starts. Somewhat 
similar behavior of the gauge parameters has appeared in the study of the exotic H^^-ion 



In Figs. 9-10 the behaviors of Req, Rmax, Rcr vs inclination at B = la.u. and 10^^ G are 
displayed. The behavior of Reg vs 9 demonstrates almost no dependence on 9 in contrast 
to both Rmax and Rcr which drastically decrease with the growth of 9. We do not have a 
reliable physical explanation of this behavior. 

The total energy dependence of H2 (at R = Reg) as a function of the inclination angle 9 
for different magnetic fields is shown in Fig. 11. The dotted line corresponds to the if-atom 
total energy in the corresponding magnetic field. For weak magnetic fields the hydrogen atom 



i 



16 




B (Gauss) 



FIG. 6: i7^-ion: domains of existence ^ non-existence for the Ig state. 

total energy is always higher than that of the if2'^-ion. However, for i? > 2 x 10^^ G the 
situation changes - a minimum of the total energy for angles 6 > 6cr does not exist any 
more. Surprisingly, 6cr corresponds approximately to the moment when the total energy of 
the if-atom becomes equal to the total energy of the if2'^-ion. If the form of vector potential 
(0) is kept fixed with ^ = 1/2 and Y = Z = 0{d = 0), then a spurious minimum appears; 
its position is displayed by the dotted curve. However, if the gauge center parameters are 
released this minimum disappears (see the discussion above). It was an underlying reason 
for the erroneous statement about the existence of the unstable H2 ion in this domain with 
a possibility to dissociate H2 H + p (see [0]). For all magnetic fields studied the total 
energy is minimal at ^ = 0° (parallel configuration) and then increases monotonically with 
inclination in complete agreement with statements of other authors |]^, [TT|, |T^, |TB[ . 

In a similar way the binding energy Eb = B — Et, as well as the dissociation energy 
(affinity to a hydrogen atom) Ed = Eh — Et as a function of 6 always decreases when 
changing from the parallel to the perpendicular configuration (see Fig. 11). Such behavior 
holds for all values of the magnetic field strength studied. Thus, we can draw the conclusion 
that the molecular ion becomes less and less stable monotonically as a function of inclination 



angle. This confirms the statement made in H, 11, 14, 16|, that the highest molecular stability 
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FIG. 7a: The dependence d vs R at B = 1 a.u. for different inclinations 9 ^ 0°. 

of the Ig state of H2 occurs for the parallel configuration. Thus, the H2 molecular ion is 
the most stable in parallel configuration. 

We extend the vahdity of this statement to magnetic field strengths lO-*^^ < B < 4.414 x 
10^^ G. It is worth emphasizing that the rate of increase of binding energy E^, with magnetic 
field growth depends on the inclination - it slows down with increased inclination. This 
effect implies that the i7^-ion in the parallel configuration becomes more and more stable 
against rotations - the energy of the lowest rotational state increases rapidly with magnetic 
field (see Table V below and the discussion there) . 

Regarding the interproton equilibrium distance R^q, one would naively expect that it 
would always decrease with inclination (see Fig. 12). Indeed, for all the magnetic fields 
studied we observe that R^q at ^ = 0° is larger than for any 6' 7^ 0° (see below. Tabs. I, II, 
III). This can be explained as a natural consequence of the much more drastic shrinking of 
the electronic cloud in the direction transverse to the magnetic field than in the longitudinal 
one. Actually, for magnetic fields B < 10^^ G the equihbrium distance Req decreases mono- 
tonically with inclination growth until it reaches 9cr, as seen in Fig. 12. As mentioned above. 
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FIG. 7b: The dependence ^ vs R at B = la.u. for different inclinations 9 ^ 0°. 



if the parameters of the vector potential are kept fixed, ^ = 1/2 and y = Z = 0((i = 0),a 
spurious minimum appears and generates anomalous (spurious) R^q behavior for 9 > 9cr{see 



In Tabs. I, II and III the numerical results for the total energy Ex, binding energy Ej, 
and equilibrium distance Reg are displayed for 9=0°, 45° and 90°, respectively. As seen 
in Table I, our results for 9 = 0° lead to the largest binding energies for B > 10^^ G in 
comparison with previous calculations. For B < 10^^ G, our binding energies for the parallel 
configuration appear to be very close (of the order of < 10""^ — 10~^ in relative deviation) 
to the variational results by Wille ||14|, which are the most accurate so far in this region of 
magnetic field strengths. Those results are based on the use of a trial function in the form 
of a linear superposition of ~ 500 Hylleraas type functions. It is quite striking that our 
simple trial function (8) with ten variational parameters gives comparable (for B < 10^^ G) 
or even better (for B > 10^^ G) accuracy. It is important to reveal the reason why the trial 
function fails to be increasingly inaccurate with magnetic field growth for B > 10^^ G. 



An explanation of this inaccuracy is related to the fact that in the {x, y)- directions the exact 
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FIG. 8a: The dependence d vs R at B = 10 G for different inclinations 9 ^ 0°. 



wave function decays asymptotically as a Gaussian function, unlike the Hylleraas functions 
which decay as the exponential of a linear function. The potential corresponding to the 
function [|1^] reproduces correctly the original potential near Coulomb singularities but fails 
to reproduce ~ p^-growth at large distances. This implies a zero radius of convergence of 
the perturbation theory for which the variational energy represents the first two terms (see 
discussion in [pS]). 

The results for 6 = 45° are shown in Table II, where a gradual shortening of the equi- 
librium distance is accompanied by an increase of total and binding energies with magnetic 
field. It is worth noting that the parameter C, evolves from about 0.5 to 0.93 with magnetic 
field growth, thus changing from the symmetric gauge for weak fields to an almost asym- 
metric one for strong ones. This phenomenon takes place for all orientations < 6 < 6cr, 
becoming more and more pronounced with increasing inclination angle (see below). We are 
unaware of any other calculations for 6 = 45° to compare ours with. 

For the perpendicular configuration 6 = 90°, the results are presented in Table III. Similar 
to what appeared for the parallel configuration (see above) our results are again slightly less 
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FIG. 8b: The dependence ^ vs R at B = 10 G for different inclinations 9 ^ 0°. 



accurate than those of Wille fl^ for B < 10^" G, but becoming the most accurate results 
for stronger fields. In particular, it indicates that the domain of applicability of a trial 
function in the form of a superposition of Hylleraas type functions, becomes smaller as the 
inclination grows. The results reported by Larsen |Tl]] and by Kappes-Schmelcher |T6| are 



slightly worse than ours, although the difference is very small. The evolution of the gauge 
parameters follow a similar trend, as was observed at ^ = 45°. In particular, ^ varies from 
^ = 0.64 to ^ = 0.98 with magnetic field growth from 5 = 10^ G to ~ 2 x 10" G [|1|]. We 
should emphasize that the results of Larsen and Wille ||T^ for B > 10^^ G do not seem 
relevant because of loss of accuracy, since the if^ ion does not exist in this region. 

In order to characterize the electronic distribution of H2 for different orientations we 
have calculated the expectation values of the transverse < p > and longitudinal < \z\ > 
sizes of the electronic cloud (see Table IV). Their ratio is always limited, 

^ < 1- 

< \z\ > 

and quickly decreases with magnetic field growth, especially for small inclination angles. 
This reflects the fact that the electronic cloud has a more and more pronounced needle-like 
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FIG. 9: The dependence of Rcrit and position of the maximum Rmax compared to the equilibrium 
position Req at B = la.u. for different inclinations 6. 



form oriented along the magnetic line, as was predicted in the classical papers [^H^l- The 
behavior of < p > itself does not display any unusual properties, smoothly decreasing with 
magnetic field, quickly approaching the cyclotron radius for small inclinations and large 
magnetic fields. In turn, the < |z| > monotonically decreases with inclination growth. 

As already mentioned, the results of our analysis of the parallel configuration of H2 
turned out to be optimal for all magnetic fields studied, being characterized by the smallest 
total energy. Therefore, it makes sense to study the lowest vibrational and also the lowest 
rotational state (see Table V). In order to do this we separate the nuclear motion along 
the molecular axis near equilibrium in the parallel configuration (vibrational motion) and 
deviation in 6 of the molecular axis from ^ = 0° (rotational motion). The vicinity of the 
minimum of the potential surface E{6, R) at 6 = 0°, R = Reg is approximated by a quadratic 
potential, and hence we arrive at a two-dimensional harmonic oscillator problem in the {R, 6)- 
plane. Corresponding curvatures near the minimum define the vibrational and rotational 
energies (for precise definitions and discussion see, for example, [|ll|])- We did not carry out a 
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FIG. 10: The dependence of Rcrit and position of the maximum Rmax compared to the equihbrium 
position R(.q at B = 10^^ G for different inchnations 9. 



detailed numerical analysis, making only rough estimates of the order of 20%. For example, 
at B= 10^^ G we obtain E^if, = 0.276 Ry in comparison with E^if, = 0.259 Ry given in [0, 
where a detailed variational analysis of the potential electronic curves was performed. Our 
estimates for the energy, -E^b, of the lowest vibrational state are in reasonable agreement 
with previous studies. In particular, we confirm a general trend of the considerable increase 



of vibrational frequency with the growth of B indicated for the first time by Larsen |]TT 



The dependence of the energy on the magnetic field is much more pronounced for the lowest 
rotational state - it grows much faster than the vibrational one with magnetic field increase. 
This implies that the ifg'^-ion in the parallel configuration becomes more stable for larger 
magnetic fields (see the discussion above). From a quantitative point of view the results 
obtained by different authors are not in good agreement. It is worth mentioning that our 



results agree for large magnetic fields > 10 a.u. with results by Le Guillou et al. [|l5l, obtained 
in the framework of the so called 'improved static approximation', but deviate drastically at 



B = 1 a.u. , being quite close to the results of Larsen |11] and Wille |13|. As for the energy of 
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FIG. 11: H2 total energy {Et) for the ground state \g as function of the inclination angle 9 for 



different magnetic fields. The dotted hues correspond to the H-atom total energy taken from |2C]. 
Dashed lines describe a total energy corresponding to a spurious minimum (see discussion in the 
text) . 
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FIG. 12: E.2 equilibrium distance as function of the inclination angle B for the 1^ state. Dashed 
lines describe the position of a spurious minimum (see discussion in the text and Fig. 11). 
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TABLE I: Total, Et, binding, 
parallel configuration, 6 = 0°. ^ 



Eh, energies and equilibrium distance R, 
This value is taken from 12 



eq 



for the state Ig in the 



B 



Et (Ry) 



Eb (Ry) 



Req (a.U.) 



B = 



10^ G 



1 a.U. 



10^0 G 



10 a.U. 



10" G 



100 a.U. 



1012 G 



1000 a.U 



10" G 



4.414 X IQi^G 



-1.20525 
-1.20527 
-1.15070 
-1.15072 
-0.94991 

-0.94642 
1.09044 
1.09031 
5.65024 



35.0434 
35.0428 

89.7090 



408.3894 

408.566 
977.2219 



4219.565 
4231.82 

18728.48 



1.57623 
1.57625 
1.94991 

1.9498 
1.94642 
3.16488 
3.16502 
4.34976 
4.35 
4.35 

4.3346 
7.50975 

7.5104 
7.34559 
10.2904 
10.2892 
10.1577 

10.270 
10.2778 
17.1425 
17.0588 

16.966 
22.7781 
21.6688 
22.7069 

22.67 
22.7694 
35.7539 

23.52 

35.74 
54.4992 



1.9971 
1.997 
1.924 
1.924 
1.752 
1.752 

1.76 
1.246 
1.246 
0.957 
0.950 
0.958 
0.950 
0.593 
0.593 

0.61 
0.448 
0.446 
0.455 
0.448 
0.446 
0.283 

0.28 
0.278 
0.220 
0.219 
0.221 
0.222 
0.219 
0.147 
0.125 

0.15 
0.101 
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the lowest rotational state, our results are in good agreement with those obtained by other 
authors (see Table V). 

In Fig. 13 we show the electronic distributions J dy\ilj{x,y, z)\'^ for magnetic fields 
10^, 10^°, 10^^, 10^^ G and different orientations for H2 in the equilibrium configuration, 
R = Reg. It was already found explicitly p2| that a.t 6 = 0° with magnetic field in- 
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TABLE II: Total E^, binding Eh energies and equilibrium distance R^q for the Ig state a,t 6 = 45°. 
Optimal value of the gauge parameter ^ is given and d = is assumed (see text) . 



B 


Et (Ry) 


Eb (Ry) 


Rpf, (a.u.) 




10^ G 


-1.14248 


1.56801 


1.891 


0.5806 


1 a.u. 


-0.918494 


1.918494 


1.667 


0.5855 


IQio G 


1.26195 


2.99337 


1.103 


0.5958 


10 a.u. 


6.02330 


3.97670 


0.812 


0.6044 


10^^ G 


36.15633 


6.39686 


0.466 


0.6252 


100 a.u. 


91.70480 


8.29520 


0.337 


0.6424 


10^2 G 


413.2987 


12.2332 


0.198 


0.6890 


1000 a.u. 


985.1956 


14.8044 


0.147 


0.7151 



FIG. 13: Electronic distributions / dy\i/}{x,y, z)\'^ (normalized to one) for the Ig state of H2 
(equilibrium configuration) for different magnetic fields and inclinations. 

crease there is a change from 'ionic' (two-peak electronic distribution) to 'covalent' cou- 
pling (one-peak distribution). We find that a similar phenomenon holds for all inclina- 
tions. If for B — 10^ G, all electronic distributions are characterized by two peaks for all 
inclinations, then for B — 10^^ G all distributions have a single sharp peak. The 'sharp- 
ness' of the peak grows with magnetic field. Fig. 13 also demonstrates how the change of 
the type of coupling appears for different inclinations ~ for larger inclinations a transition 
{two — peaks) (one — peak) appears for smaller magnetic fields. It seems natural that 
for the perpendicular configuration 9 = 90°, where the equilibrium distance is the smallest, 
this change appears for even smaller magnetic field. 

In Figs. 14-18 we present the evolution of the electronic distributions as a function of 
interproton distance R, for inchnations 0°, 45° a,t B — la.u. and 10^^ G together with the 
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TABLE III: Total, Ex, and binding, E^, energies and the equilibrium distance R^q for the Ig 
state in the perpendicular configuration, 9 = 90°. The optimal value of the gauge parameter ^ is 
presented and d is kept fixed, d = (see text). 



B 



Et (Ry) 



Eb (Ry) 



Req (a.U. 



10^ G 



1 a.U. 



W G 



10 a.U. 



lOii G 



100 a.U. 



1012 G 



1000 a.U. 



-1.137342 

-0.89911 

-0.89774 
1.36207 

6.23170 

36.7687 

92.7346 



1.56287 
1.56384 
1.89911 
1.8988 

I. 8977 
2.89324 
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3.76830 
3.7620 
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5.6818 
7.26543 

7.229 

4.558 

II. 58 



1.875 
1.879 
1.635 
1.634 
1.65 
1.059 
1.067 
0.772 
0.772 
0.442 
0.428 
0.320 
0.320 

0.148 

0.1578 
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0.7329 
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i?-dependence for the inclination 90° aX B = la.u. . The values of magnetic fields are chosen 
to illustrate in the most explicit way the situation. In all figures a similar picture is seen. 
Namely, at not very large magnetic fields B < 10^^ G and for all inclinations the electronic 
distribution at small R < Rcr is permutationally symmetric and evolves with increase of R 
from a one-peak to a two-peak picture with more and more clearly pronounced, separated 
peaks. Then for R = Rcr this symmetry becomes broken and the electron randomly chooses 
one of protons and prefers to stay in its vicinity. For R ^ R^r the electronic distribution 
becomes totally asymmetric, the electron looses its memory of the second proton. This 
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TABLE IV: The Ig state: Expectation values of the transverse < p > and longitudinal < |2;| > 
sizes of the electron distribution for the H^-ion in a.u. at different orientations and magnetic field 
strengths. At ^ = 0° the expectation value < p > almost coincides to the cyclotron radius of the 
electron. 



B 


0" 


<P> 
45° 


90° 


0° 


< k > 
45° 


90° 


10^ G 


0.909 


1.002 


1.084 


1.666 


1.440 


1.180 


1 a.u. 


0.801 


0.866 


0.929 


1.534 


1.313 


1.090 


IQio G 


0.511 


0.538 


0.569 


1.144 


0.972 


0.848 


10 a.u. 
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IQii G 
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0.009 
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signals that the chemical reaction H2 ^ H + p has already happened. For larger magnetic 
fields B > IQii G for i? < Rcr the electronic distribution is always single-peaked, a transition 
from a one-peak to a two-peak picture occurs for R > R^r, where the electronic distribution 
is already asymmetric. 



FIG. 14: Evolution of the electronic distributions J dy\'tp{x,y, z)\'^ (normalized to one) and their 
contours for the Ig state of the (ppe) system with interproton distance for B = 1 a.u., 9 = 0°. 

To complete the study of the Ig state we show in Fig. 19 the behavior of the variational 
parameters of the trial function (|^) as a function of the magnetic field strength for the 
optimal (parallel) configuration, 6 = 0°. In general, the behavior of the parameters is rather 
smooth and very slowly-changing, even though the magnetic field changes by several orders 
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TABLE V: Energies of the lowest vibrational {Emb) and rotational (Erot) electronic states associ- 
ated with the Ig state at ^ = 0°. The indexes in Le Guillou et al |jl^ are assigned to the 'improved 
adiabatic approximation' (a), and to the 'improved static approximation' (b). 

B Et (Ry) E,,t (Ry) Erot (Ry) 



109 G 



1 a.u. 



IQio G 



10 a.u. 



10" G 
100 a.u. 
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35.0434 
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FIG. 15: Evolution of the electronic distributions / dy\ip{x,y, z)\'^ (normalized to one) and their 
contours for the Ig state of the (ppe) system with interproton distance for B = 1 a.u., 9 = 45°. 
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FIG. 16: Evolution of the electronic distributions f dy\tp{x,y, z)\'^ (normalized to one) and their 
contours for the Ig state of the (ppe) system with interproton distance for B = 1 a.u., 6 = 90°. 



FIG. 17: Evolution of the electronic distributions / dy\ip{x,y, z)\'^ (normalized to one) and their 
contours for the Ig state of the (ppe) system with interproton distance for B = 10^^ G,9 = 0°. 

of magnitude. This is in drastic contrast with the results of Kappes-Schmelcher (see 
Fig. 1 in this paper). In our opinion such behavior of the parameters of our trial function (9) 
reflects the level of adequacy (or, in other words, indicates the quahty) of the trial function. 
In practice, the parameters can be approximated by the spline method and then can be used 
to study magnetic field strengths other than those presented here. 



FIG. 18: Evolution of the electronic distributions / dy\ip{x,y, z)\'^ (normalized to one) and their 
contours for the Ig state of the (ppe) system with interproton distance for B = 10^^ G,9 = 45°. 
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FIG. 19: Variational parameters of the trial function (9) as a function of the magnetic field strength 



B for the Ig state in the parallel configuration, 9 = 0°. In this case /?i = 



Plx _ Ply [3 _ §2x 
2 5 Mi 2 



^ - ^ (see eqs. (5,6,8), cf. 



'3 ~ ^ ~ ~2^ ^'^^^ eqs.i^o,D,o;, ci. [|£^j The parameter A3 is fixed, being equal to 1 and 
^ = 1/2, (i = (see text). The error bars correspond to relative deviation in the variational energy 
in the region AE^t = < 10"^ 
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V. CONCLUSION 



We have carried out an accurate, non-relativistic calculation in the Born-Oppenheimer 
approximation for the lowest state of the H2 molecular ion for different orientations of the 
magnetic field with respect to the molecular axis. We studied constant uniform magnetic 
fields ranging from zero up to -B = 4.414 x 10^'^ G, where non-relativistic consideration holds. 

For all magnetic fields studied there exist a region of inclinations for which a well- 
pronounced minimum in the total energy surface for the Ig state of the system {ppe) is found. 
This shows the existence of the H2 molecular ion for magnetic fields 5 = — 4.414 x 10^^ G. 
The smallest total energy is always found to correspond to the parallel configuration, ^ = 0°, 
where protons are situated along the magnetic line. The total energy increases, while the 
binding energy decreases monotonically as the inclination angle grows. The rate of total 
energy increase as well as binding energy decrease is seen to be always maximal for the par- 
allel configuration. The equilibrium distance exhibits quite natural behavior as a function 
of the orientation angle 9 - for fixed magnetic field the shorter equilibrium distance always 
corresponds to the larger 9. 

Confirming the qualitative observations made by Khersonskij 0] for the Ig state in the 
contrast to statements in [[ll| , , we demonstrate accurately that the H2-ion does not exist 
at a certain range of orientations for magnetic fields B >2 x 10^^ G. As the magnetic field 
increases the region of inclinations where H2 does not exist is seen to broaden, reaching 
rather large domain 25° < 9 ^ 90° for B = 4AU x 10^^ G. 

We find that the electronic distributions for H2 in the equilibrium position are qualita- 
tively different for weak and large magnetic fields. In the domain B < 10^^ G the electronic 
distribution for any inclination has a two-peak form, peaking near the position of each pro- 
ton. On the contrary for B > 10^^ G the electronic distribution always has a single peak 
form with the peak near the midpoint between the protons for any inclination. This implies 
a physically different structure for the ground state - for weak fields the ground state can 
be modelled as a 'superposition' of hydrogen atom and proton, while for strong fields such 
modelling is not appropriate. 

Unlike standard potential curves for molecular systems in the field-free case, we observe 
for > 0" that each curve has a maximum and approaches to the asymptotics at R 00 
from above. The electronic distribution evolves with R from a one-peak form at small R 
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to a two-peak one at large R. There exists a certain critical Rcr at which one of peaks 
starts to diminish, manifesting a breaking of permutation symmetry between the protons 
and simultaneously the beginning of the chemical reaction H2 H + p. 

Combining all the above-mentioned observations we conclude that for magnetic fields of 
the order of magnitude B ~ 10^^ G some qualitative changes in the behavior of the H2 ion 
take place. The behavior of the variational parameters also favors this conclusion. This 
hints at the appearance of a new scale in the problem. It might be interpreted as a signal 



of a transition to the domain of developed quantum chaos (see, for example, |33]). 
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